Bal = input(' Balanced Load? (1=Yes,2=No) ');
% 用户输入负载是否平衡的选项，1代表平衡负载，2代表非平衡负载。

Dir = input(' Conversion Direction?(1=D->Y,2=Y->D)');
% 用户输入转换的方向，1代表Delta到WYE转换，2代表WYE到Delta转换。

if Dir == 1
    disp(' DELTA to WYE CONVERSION');
    % 显示Delta到WYE转换的标题。
    
    if Bal == 1
        magZab = input(' magZab= ');
        magZbc = magZab;
        magZca = magZab;
        angZab = input(' angZab= ')*pi/180;
        angZbc = angZab;
        angZca = angZab;
        % 对于平衡负载，用户输入Zab和Zbc的模值magZab、Zca和Zbc的角度angZab。
        
    elseif Bal == 2
        magZab = input(' magZab= ');
        angZab = input(' angZab= ')*pi/180;
        magZbc = input(' magZbc= ');
        angZbc = input(' angZbc= ')*pi/180;
        magZca = input(' magZca= ');
        angZca = input(' angZca= ')*pi/180;
        % 对于非平衡负载，用户输入Zab、Zbc和Zca的模值和角度。
    end
    
    Zab = magZab * exp(j * angZab);
    Zbc = magZbc * exp(j * angZbc);
    Zca = magZca * exp(j * angZbc);
    den = Zab + Zbc + Zca;
    Za = Zab * Zca / den;
    Zb = Zbc * Zab / den;
    Zc = Zca * Zbc / den;
    % 计算WYE网络中的阻抗Za、Zb和Zc。
    
    disp(' ');
    disp(' WYE IMPEDANCES in POLAR FORM');
    disp([abs(Za) angle(Za)*180/pi; abs(Zb) angle(Zb)*180/pi; abs(Zc) angle(Zc)*180/pi]);
    % 显示WYE网络中阻抗的模值和角度。
    
elseif Dir == 2
    disp(' WYE to DELTA CONVERSION');
    % 显示WYE到Delta转换的标题。
    
    if Bal == 1
        magZa = input(' magZa= ');
        magZb = magZa;
        magZc = magZa;
        angZa = input(' angZa= ')*pi/180;
        angZb = angZa;
        angZc = angZa;
        % 对于平衡负载，用户输入Za、Zb和Zc的模值magZa以及角度angZa。
        
    elseif Bal == 2
        magZa = input(' magZa= ');
        angZa = input(' angZa= ')*pi/180;
        magZb = input(' magZb= ');
        angZb = input(' angZb= ')*pi/180;
        magZc = input(' magZc= ');
        angZc = input(' angZc= ')*pi/180;
        % 对于非平衡负载，用户输入Za、Zb和Zc的模值和角度。
    end
    
    Za = magZa * exp(j * angZa);
    Zb = magZb * exp(j * angZb);
    Zc = magZc * exp(j * angZc);
    num = Za * Zb + Zb * Zc + Zc * Za;
    Zab = num / Zc;
    Zbc = num / Za;
    Zca = num / Zb;
    % 计算Delta网络中的阻抗Zab、Zbc和Zca。
    
    disp(' ');
    disp(' DELTA IMPEDANCES in POLAR FORM');
    disp([abs(Zab) angle(Zab)*180/pi; abs(Zbc) angle(Zbc)*180/pi; abs(Zca) angle(Zca)*180/pi]);
    % 显示Delta网络中阻抗的模值和角度。
end